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Abstract 

In  order  to  utilize  the  novel  electrical,  magnetic,  op¬ 
tical,  and  physical  properties  of  coated  metal  nanoparti¬ 
cles,  one  must  be  able  to  efficiently  predict  the  nanoparti¬ 
cle  size-dependent  properties  and  to  fabricate  consistently 
sized  nanoparticles.  Both  of  these  goals  can  be  obtained 
through  the  use  of  numerical  simulation.  In  this  work  the 
gold  nanoparticle  and  nanoparticle/self-assembled  mono- 
layer  systems  have  been  analyzed  with  a  physically  ac¬ 
curate  and  computationally  efficient  numerical  simulation. 
The  simulation  model  and  method  are  described  in  the  fol¬ 
lowing  paper. 


1.  Introduction 


In  this  work  we  have  used  molecular  dynamics  (MD) 
simulations  to  study  gold  nanoparticles  with  alkanethiolate 
SAMs.  This  material  system  is  an  extensively  researched 
system  in  the  study  of  SAMs  because  of  its  wide  applicabil¬ 
ity  and  the  wealth  of  experimental  data  available.  The  MD 
algorithm  is  used  because  it  can  efficiently  simulate  large 
numbers  of  atoms  with  the  desired  accuracy.  This  is  an  im¬ 
portant  consideration  when  large  numbers  of  atoms  are  to 
be  analyzed.  We  have  combined  accurate  MD  potentials 
for  each  of  the  subsystems,  namely  the  gold  nanoparticle 
and  the  alkanethiolate  polymer  chain.  The  system  is  then 
analyzed  by  varying  parameters  such  as  temperature,  SAM 
chain  length,  and  the  density  of  adsorbed  chains. 
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2.  Computational  Procedure 

Molecular  dynamics  (MD)  studies  of  alkanethiolate 
SAMs  on  gold  surfaces  are  well  documented  [9,24,29]. 
These  studies  have  been  limited  to  flat  gold  surfaces  or  fixed 
gold  atoms  in  a  nanocrystallite  [12, 15].  In  addition  to,  and 
in  support  of  these  previous  MD  simulations,  quantum  me¬ 
chanical  simulations  have  been  used  to  study  gold  clusters 
with  short  alkanethiolate  SAMs.  Interestingly,  the  quan¬ 
tum  mechanical  findings  have  suggested  that  the  gold  crys¬ 
tal  lattice  is  appreciably  affected  by  the  adsorbed  ligand  [7] 
with  an  increase  in  the  Au-Au  bond  length  of  up  to  20%. 
The  quantum  mechanical  results  suggest  that  the  MD  sim¬ 
ulations  should  include  dynamic  gold  atoms  in  order  to  ac¬ 
curately  model  the  nanoparticle/SAM  system.  Without  the 
simulation  of  dynamic  gold  atoms,  the  changes  that  occur 
to  the  surface  is  not  considered. 

2.1.  Software  and  Optimization 

The  software  for  this  project  was  chosen  based  on  its 
current  capabilities,  extensibility,  and  parallel  performance. 
The  software  used  needs  to  be  extensible  because  of  the 
many  unique  computational  algorithms  required  in  this 
project.  Some  of  the  unique  algorithms  include  comput¬ 
ing  the  radial  pressure  distribution  inside  of  the  nanoparti¬ 
cle  and  use  of  uncommon  potentials.  For  MD  simulations 
it  was  determined  that  the  Large-scale  Atomic/Molecular 
Massively  Parallel  Simulator  (LAMMPS)  from  Sandia  Na¬ 
tional  Laboratories  was  the  best  choice  as  it  is  open  source, 
widely  used,  and  has  massively  parallel  capabilities  [19]. 
LAMMPS  is  open  source  software  allowing  for  the  easy  in¬ 
tegration  of  new  routines  to  the  application,  thus  providing 
the  required  extensibility. 
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Although  LAMMPS  is  massively  parallel,  and  easy  to 
use,  there  appeared  to  be  some  inefficiencies  in  the  soft¬ 
ware,  that  if  corrected,  could  provide  improved  perfor¬ 
mance.  A  PET  (User  Productivity  Enhancement  and  Tech¬ 
nology  Transfer)  project  at  the  U.S.  Army  Research  Lab¬ 
oratory  was  underway  at  the  time  of  this  effort  to  look  at 
this  issue.  The  project  team  discovered  many  opportunities 
for  performance  enhancements  through  extensive  profiling 
and  benchmarking.  Each  of  the  following  enhancements 
provided  some  improvement  but  together  they  provided  a 
3.5x  improvement  in  run  time  for  the  alkanethiolate  SAM 
coated  gold  nanoparticle  model.  Some  of  the  biggest  im¬ 
provements  came  from  reducing  array  accesses,  using  tem¬ 
porary  variables,  restructuring  of  if-else  blocks,  and  elim¬ 
inating  conditional  blocks  by  using  C++  templates.  Some 
other  improvements  came  from  machine  specific  changes 
to  the  code  that  assisted  the  compiler  such  as  the  manual 
unrolling  of  inner  loops  by  four. 

2.2.  Simulating  the  Gold  Substrate 

When  explicitly  considering  the  gold  substrate  as  dy¬ 
namic  atoms,  the  embedded  atom  method  (EAM)  [4, 5, 18] 
and  the  effective  medium  (EM)  [11,20,25]  potentials  have 
received  the  most  attention,  are  computationally  efficient, 
and  provide  accurate  results  in  MD  simulations  under  many 
conditions.  These  methods  use  empirical  embedding  func¬ 
tions  that  originate  from  density  functional  theory  (DFT) 
and  have  been  successfully  used  to  model  many  of  the  phys¬ 
ical  properties  of  gold.  Since  the  embedding  functions  are 
empirical,  they  are  computationally  efficient  and  yet  still  ac¬ 
curate  enough  for  many  purposes.  This  efficiency  makes  it 
possible  for  large  numbers  of  atoms  to  be  modeled,  and  thus 
analysis  of  realistically  sized  systems  is  feasible.  We  have 
decided  to  use  the  EAM  for  modeling  the  gold  substrate  be¬ 
cause  of  its  accuracy  and  the  extensions  available  for  mod¬ 
eling  surfaces,  such  as  the  extended  EAM  (XEAM)  [13]. 
In  the  EAM,  the  potential  function  is  a  combination  of  the 
embedding  energy  and  pair  potential  terms  as  in  equation  1 . 
The  actual  data  that  is  used  in  this  work  is  given  in  Figure 
1  along  with  a  second  dataset  from  Voter  and  Chen  [28].  A 
detailed  discussion  of  the  method  as  used  in  this  work  can 
be  found  in  Daw  et  al.  [4,5]. 

E tot  -  /*  Fi  ( Ph,i )  +  /*  '  'y  '  <$>ij  (Rij)  (1) 

*  *  i(^i) 

The  gold  nanoparticle  model  begins  as  a  simple  sphere 
with  a  face  centered  cubic  (FCC)  lattice  structure  as  seen  in 
figure  2(a).  The  model  is  then  heated  to  1000  K,  or  above 
the  melting  temperature,  and  held  for  a  short  time  of  j  10 
ns.  After  the  nanoparticle  has  been  melted  for  about  10  ns 


(a)  Electron  density 


(b)  Embedding  energy  where  pi  has  been 
scaled  by  p &  =  0.025  for  auu3  and  0.15  for 
Voter-Chen. 


(c)  Pair  potential 

Figure  1.  Gold  (Au)  EAM  potentials  from 
LAMMPS  [19].  Distances  are  measured  in 
Angstroms  (A)  and  energy  is  given  in  elec¬ 
tron  volts  (eV) 
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(a)  Unequilibrated  (b)  Equilibrated 


Figure  2.  Equilibrated  and  Unequilibrated 

the  model  is  slowly  cooled  to  the  desired  simulation  tem¬ 
perature.  This  process  minimizes  the  energy  of  the  sys¬ 
tem  and  provides  the  most  accurate  representation  of  the 
equilibrated  nanoparticle  possible.  The  equilibrated  gold 
nanoparticle  in  figure  2(b)  has  large  facets  and  is  no  longer 
spherical.  This  geometry  is  expected  as  it  minimizes  the 
energy  of  the  system  as  opposed  to  the  perfect  sphere. 

2.3.  Modeling  the  Alkanethiolate  Chain 

SAMs  are  formed  by  the  spontaneous  adsorption  of  self¬ 
assembling  molecules  on  a  solid  substrate.  A  SAM’s  ther¬ 
mal  stability  and  structural  properties  are  governed  by  the 
substrate-head  molecular  and  intrachain  interactions  in  the 
chemisorbed  layer.  The  chemical  identity  of  the  SAM  tail 
molecule  governs  the  surface  properties  of  the  coated  metal 
and  affects  interactions  with  the  environment.  Examples 
of  SAM  systems  include  alkanethiols  or  dialkyl  disulfides 
on  gold  and  silver,  organosilicon  on  hydroxylated  silicon 
wafers,  and  carboxylic  acids  on  aluminum. 

The  alkanethiolate  polymer  chain  that  coats  the  gold 
nanoparticle  in  this  work  may  be  modeled  with  a  range  of 
complexity.  For  MD  simulations  the  three  methods  avail¬ 
able  are  the  high-accuracy/high  computational  requirement 
all-atom  model,  the  less  accurate/lower  computational  re¬ 
quirement  united  atom  method,  and  the  coarse  grained  bond 
model.  For  instance,  consider  the  CH2  molecule  that  makes 
up  the  backbone  of  the  alkanethiolate  chain.  When  us¬ 
ing  the  all-atom  model  the  numerical  simulation  must  con¬ 
sider  three  particles  for  each  molecule.  Using  the  united 
atom  model  only  one  particle  is  simulated  for  each  molecule 
[10],  and  with  the  coarse  grained  model  multiple  molecules 
are  grouped  together  into  one  particle  with  perhaps  five 
or  more  CH2  molecules  grouped  into  one  simulated  par¬ 
ticle  [6,  17,27],  The  model  chosen  for  a  particular  simu¬ 
lation  depends  greatly  upon  the  level  of  detail  required  in 
the  simulation  and  corresponding  results  of  interest.  For  the 


Figure  3.  Sample  alkanethiolate  chain  with  8 
carbon  atoms  along  the  backbone. 


simulations  in  this  work  we  have  chosen  the  united  atom 
method,  which  is  computationally  more  efficient  that  the  all 
atom  method  yet  provides  enough  accuracy  for  any  antici¬ 
pated  quantitative  analysis. 

The  modeling  of  the  alkanethiolate  chain  using  the 
united  atom  model  is  consistent  with  many  of  the  Au-SAM 
system  published  results.  An  example  of  the  alkanethiolate 
polymer  chain  model  used  in  this  work  is  shown  in  Figure 
3.  The  potentials  used  in  the  united  atom  model  include  a 
bending  potential  for  the  S-C-C  bond  and  the  C-C-C  bond. 
The  bending  potential  can  be  modeled  as  either  a  flexible 
bond  with  a  form  similar  to  equation  2  or  as  a  stiff  bond 
where  the  angle  is  held  constant.  There  is  a  dihedral  poten¬ 
tial  for  X-C-C-X  bonds,  where  X  can  be  either  S  or  C  and 
takes  the  form  of  equation  3.  There  are  two  common  meth¬ 
ods  used  in  the  literature  to  model  the  bond  lengths  along 
the  chain.  The  bond  lengths  along  the  chain  can  be  held 
constant  [24]  with  the  RATTLE  [2]  algorithm  or  they  can 
be  modeled  as  harmonic  bonds  [17,21]  as  given  in  equa¬ 
tion  4.  Using  a  harmonic  bond  to  describe  the  C-C  and  S-C 
backbone  bonds  allows  the  distance  between  the  backbone 
atoms  to  change  while  not  allowing  the  bonds  to  break. 


Ubending  (#)  =  ^kg  (9  -  90)2 

(2) 

5 

U dihedral  (0)  —  ^  ^  (XiCOS  (</>) 

(3) 

2  —  0 

Ubond  ( r )  =  \kr  (r  -  r0f 

(4) 

The  values  of  the  parameters  used  in  this  work  are  as 
follows.  For  the  C-C  bond  stretch  model  kr  =  30.3492]-, 

ro  =  1.53A,  and  for  S-C  bond  stretching  kr  =  35.594-2]-, 

A 

ro  =  1 .82,4  [21],  The  bond  bending  potential  for  the  C- 
C-C  angle  is  kg  =  8840.0419^,  90  =  109.5  deg,  and 
for  the  S-C-C  angle  kg  =  8840.0419^,  90  =  114.4  deg 
[24],  The  parameters  for  the  dihedral  bonds  from  equation 
3  are  a0  =  0.0961661,  oq  =  0.125981,  a2  =  -0.1359773, 
a3  =  -0.0317102,  a4  =  0.271954,  a5  =  -0.326414  [24], 
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2.4.  Binding  of  Sulfur  to  Gold 


The  binding  of  the  head  sulfur  atom  of  an  alkanethiolate 
chain  to  the  gold  substrate  is  of  great  interest  as  it  affects  the 
location,  orientation,  movement,  and  desorption  of  chains 
from  the  substrate  [3,7, 16],  Much  effort  has  been  devoted  to 
finding  accurate  potentials  for  simulating  the  sulfur  binding 
energy.  Determining  the  potential  parameters  in  MD  sim¬ 
ulations  requires  the  use  of  expensive  ab  initio  calculations 
such  as  DFT  [3,7].  The  most  published  potentials  resulting 
from  these  calculations  are  the  12-3  potential  [8],  used  most 
commonly  in  Monte  Carlo  (MC)  simulations  of  alkanethio- 
lates  on  gold  [9, 24]  and  the  Morse  potential,  used  in  many 
of  the  MD  simulation  studies  [9,  14,  16,29].  In  addition 
to  the  issue  of  which  potential  best  describes  the  binding 
between  SAM  head-group  and  substrate,  the  possibility  of 
SAM  mobility  is  also  an  important  consideration  [3],  On 
curved  surfaces,  as  opposed  to  flat  surfaces,  the  boundary 
conditions  that  restrict  lateral  movement  on  the  surface  are 
not  present  and  so  chain  mobility  is  a  primary  consideration. 

3.  Computational  Results 

In  this  section  we  will  present  our  findings  from  the  MD 
simulations  studied.  Of  interest,  is  the  observation  that  for 
gold  nanoparticles,  adsorbed  SAMs  modify  the  crystallinity 
of  the  gold  atoms  near  the  nanoparticle  surface.  First,  the 
structural  change  needs  to  be  quantified,  and  then  the  abil¬ 
ity  to  predict  the  structural  change  from  parameters  such  as 
temperature,  alkanethiolate  chain  length,  surface  coverage, 
etc.  can  be  investigated.  The  following  discussions  will 
focus  on  the  observed  properties  of  the  simulated  nanopar¬ 
ticle. 

In  order  to  study  the  affects  of  the  SAM  on  the  nanopar¬ 
ticle  we  have  computed  the  pair  correlation  function,  diffu¬ 
sion  coefficients,  and  radial  pressure  and  density  distribu¬ 
tions.  Each  of  these  results  will  give  some  insight  into  what 
changes  the  nanoparticle  surface  experiences  by  the  addi¬ 
tion  of  the  adsorbed  SAMs.  For  instance,  the  diffusion  co¬ 
efficient  will  be  used  to  compute  the  mobility  of  the  SAMs 
that  are  adsorbed  on  the  surface,  and  also  the  mobility  of 
the  gold  atoms  near  the  surface  of  the  nanoparticle.  This  in¬ 
formation  will  indicate  the  phase  of  the  materials  at  and  on 
the  surface  and  also  help  in  predicting  if  nucleation  of  the 
SAMs  may  occur,  where  under  low  surface  coverage  condi¬ 
tions  the  chains  will  form  dense  groups  and  uneven  surface 
coverage.  The  density  distribution  is  used  to  find  the  depth 
to  which  the  sulfur  atoms  penetrate  the  surface  of  the  gold 
nanoparticle.  The  density  distribution  is  used  to  define  a 
corrugation  factor  that  describes  the  change  in  the  size  of 
the  nanoparticle  core. 


Figure  4.  Pair  correlation  function  for  bare 
5nm  gold  nanoparticle.  Note  the  change  in 
shape  and  number  of  peaks  between  900  K 
and  1000  K.  This  change  indicates  that  a 
phase  change  from  solid  to  liquid  has  taken 
place. 


3.1.  Pair  Correlation  Function 

The  pair  correlation  function  g  (r)  is  the  number  of 
atoms  a  distance  r  from  a  given  atom  compared  with  the 
number  at  the  same  distance  in  an  ideal  gas  at  the  same  den¬ 
sity  [1].  The  pair  correlation  therefore  is  used  to  provide  in¬ 
sight  into  the  structure  of  the  nanoparticle,  from  highly  crys¬ 
talline  with  sharp  peaks  for  nearest  neighbors,  next-nearest 
neighbors,  and  so  on,  to  liquids  with  peaks  at  multiples  of 
the  nearest  neighbor  distance.  Equation  5  is  used  to  com¬ 
pute  the  pair  correlation  function  g  (r). 

=  (5) 

\  i  jj^i  / 

In  equation  5  V  is  the  volume  of  the  simulation  bounds, 
N  is  the  number  of  simulated  particles,  and  6  is  the  delta 
function. 

The  pair  correlation  function  is  useful  because  the  en¬ 
semble  average  of  any  pair  function  may  be  expressed  as  a 
function  of  the  pair  correlation  function.  Some  results  for 
the  bare  5nm  gold  nanoparticle  are  given  in  figure  4.  These 
results  indicate  a  phase  transition  occurs  in  the  nanoparticle 
from  a  crystalline  solid  at  900  K  to  a  liquid  at  1000  K.  It  can 
therefore  be  assumed  that  a  5nm  gold  nanoparticle  melts  in 
this  temperature  range. 
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3.2.  Diffusion 


The  diffusion  coefficients  computed  here  are  used  to 
quantify  the  mobility  of  the  SAMs  on  the  gold  nanoparti¬ 
cle  surface.  The  diffusion  coefficient  of  the  alkanethiolate 
chain  is  computed  by  considering  the  mean  squared  dis¬ 
tance  traveled  by  the  sulfur  head  group  during  a  simulation. 
The  diffusion  coefficient,  D,  is  computed  using  equation  6. 


d(r2  ( t )) 
dt 


2  dD 


(6) 


In  equation  6,  d  is  the  number  of  dimensions  available 
for  atomic  diffusion,  in  this  case  3,  t  is  elapsed  time,  and 
<r2  (f)>  is  the  mean  squared  displacement  (MSD)  of  the 
atoms  being  tracked. 

In  addition  to  studying  the  diffusion  of  the  sulfur  atoms 
we  have  also  considered  the  diffusion  of  the  gold  surface 
atoms  that  are  bound  to  the  sulfur  atoms.  This  information 
is  used  to  determine  whether  or  not  the  sulfur  atoms  are 
moving  over  the  surface  atoms  or  with  the  surface  atoms. 
If  the  sulfur  is  moving  over  the  nanoparticle  surface,  they 
are  considered  mobile  and  may  nucleate  on  the  surface.  If 
the  gold  surface  is  diffusing  at  a  similar  rate  then  the  gold 
on  the  surface  may  be  carrying  the  sulfur  with  it,  and  the 
chains  are  not  moving  in  relation  to  the  gold  surface.  Fig¬ 
ure  5  contains  the  MSD  results  for  alkanethiolate  SAMs  ad¬ 
sorbed  on  a  gold  nanoparticle  at  300  K  and  full  surface  cov¬ 
erage.  These  results  give  a  MSD  value  for  the  sulfur  head 
group  of  8.05-10~7c?n/s2  and  for  the  gold  surface  atoms  of 
9.75  •  10~8c?n/s2.  The  difference  in  diffusion  coefficients 
is  about  one  order  of  magnitude,  indicating  that  the  alka¬ 
nethiolate  chains  are  relatively  mobile  over  the  nanoparticle 
surface. 

Another  study  performed  using  the  calculation  of  diffu¬ 
sion  coefficients  is  the  temperature  dependence  of  the  diffu- 
sivity  vs.  temperature.  By  comparing  the  logarithm  of  the 
diffusion  coefficient  of  the  SAMs  for  various  temperatures 
it  should  be  possible  to  see  regions  of  linear  increases  in 
In  (D)  vs.  TL i.  A  linear  curve  would  indicate  an  Arrhe¬ 
nius  distribution.  In  figure  6  there  are  three  distinct  linear 
regions  with  slopes  of  —260.9,  —461.6,  and  —4565.3,  indi¬ 
cating  three  phases  in  the  alkanethiolate  SAM. 


3.3.  Radial  Density  Distribution 


Figure  5.  The  MSD  of  sulfur  head  group  atoms 
and  the  gold  surface  atoms  in  a  SAM  coated 
nanoparticle.  These  results  indicate  that  the 
SAMs  are  moving  freely  over  the  gold  surface. 


Thus  far  the  computed  results  have  shown  the  the  5nm 
bare  nanoparticle  has  a  melting  temperature  between  900  K 
and  1000  K  (pair  correlation  function),  the  SAMs  are  mo¬ 
bile  over  the  gold  surface  (diffusion  coefficient),  and  that 
there  are  three  possible  phases  between  100  K  and  800  K. 
The  next  calculation,  namely  the  radial  density,  is  required 
as  part  of  the  radial  pressure  calculation.  The  radial  density 
distribution  is  used  to  determine  where  the  edge  of  the  gold 


Figure  6.  The  activation  energy  for  SAM  mo¬ 
bility  is  estimated  by  the  slope  of  In  (D)  vs. 
1/T  for  various  temperatures  between  100  K 
and  800  K. 
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Figure  7.  Radial  density  of  gold  in  alkanethio- 
late  SAM  coated  gold  nanoparticle  for  various 
surface  coverage  amounts. 


nanoparticle  core  occurs.  For  bare  nanoparticles  the  change 
in  density  will  be  abrupt  at  the  surface,  but  for  coated  parti¬ 
cles  the  slope  will  be  lower  if  the  SAMs  penetrate  into  the 
gold  surface.  The  SAM  penetrating  into  the  surface  was 
not  initially  assumed  because  past  work  simulating  alka- 
nethiolate  coated  gold  nanoparticles  has  assumed  that  the 
gold  atoms  do  not  move  because  of  the  interactions  with 
the  SAMs  [15].  With  the  parallel  computing  capabilities  of 
LAMMPS  it  was  possible  to  model  the  SAMs  and  the  gold 
atoms  as  dynamic  particles.  This  allowed  for  an  investiga¬ 
tion  of  whether  or  not  the  gold  atoms  are  affected  by  the 
SAMs.  When  computing  the  radial  density  as  shown  in  fig¬ 
ure  7  is  became  apparent  that  the  atoms  near  the  surface  of 
the  gold  nanoparticle  are  in  fact  moving  and  corrugating  the 
surface. 

In  an  attempt  to  quantify  the  surface  corrugation  we  have 
defined  a  corrugation  factor,  C,  that  is  defined  as 

c=  r  (P  =  0-85 Pbuik)bare  ~r(p  =  0.85 pbuik)coated 
r(P  =  0.85  Pbuik)bare 

(7) 

where  r  (p  =  0.85 pbuik)  is  the  radius  where  the  computed 
gold  density  is  85%  of  the  density  of  bulk  gold.  This  defi¬ 
nition  of  the  corrugation  factor  will  be  0.0  for  the  bare  gold 
nanoparticle  and  up  to  1 .0  for  a  gold  nanoparticle  with  im¬ 
purities  or  a  dissolved  solute.  For  coated  gold  nanoparticles 
with  various  percentages  of  surface  coverage  the  corruga¬ 
tion  factor  has  been  computed  and  is  listed  in  table  1 . 


Table  1.  Computed  corrugation  factors  for 
alkanethiolate  SAM  coated  gold  nanoparti¬ 
cles. 


Coverage 

C 

Bare 

0.000 

10% 

0.004 

20% 

0.004 

70% 

0.012 

80% 

0.018 

90% 

0.043 

100% 

0.107 

3.4.  Radial  Pressure  Distribution 

It  is  known  that  small  droplets  have  an  elevated  inter¬ 
nal  pressure  due  to  compressive  forces  arising  from  the 
surface  tension.  In  these  MD  simulations  the  radial  pres¬ 
sure  distribution  inside  of  the  nanoparticle  is  computed  us¬ 
ing  the  Irving-Kirkwood  pressure  tensor  [26].  The  Irving  - 
Kirkwood  pressure  tensor,  equation  8,  has  contributions 
from  two  terms,  namely  the  kinetic  (Pk)  and  configura¬ 
tional  (Pu). 

Pn  (r)  =  PK  (r)  +  Pu  (r)  (8) 

PN  indicates  the  normal  component  of  the  pressure  tensor. 
The  kinetic  term  is  dependent  upon  the  radial  density  distri¬ 
bution,  p  (r),  and  the  system  temperature,  equation  9. 

PK  (r)  =  kBT p  (r)  (9) 

where  ks  is  Boltzmann’s  constant.  The  configurational 
term  is  dependent  upon  atomic  interactions  that  act  across 
the  radius  of  interest,  as  seen  in  figure  8.  In  equation  10  the 
calculation  of  the  configurational  term  is  defined. 

Pu  (r)  =  -  (47rr3)  I r  •  r%|  _LduJr‘j)  (10) 

V  r«  ,lr'j 

In  equation  10  :rl3  is  the  vector  between  the  interacting 
atoms  that  crosses  radius  r  of  the  sphere  of  interest,  u  ( ) 
is  the  potential  energy  of  the  atomic  interaction  and  so 
du}r'3  ^  is  the  force  of  the  interaction.  Using  this  calculation, 
the  radial  pressure  distribution  for  bare  5nm  gold  nanopar¬ 
ticles  at  temperatures  between  200  K  and  1000  K  have  been 
computed  and  plotted  in  figure  9. 

The  pressure  results  in  figure  9  show  that  the  pressure 
at  the  surface  of  a  nanoparticle  slightly  above  the  melting 
temperature  is  about  5900  bar.  This  result  is  used  with 
the  Young -Laplace  equation,  equation  11,  to  compute  the 
nanoparticle  surface  tension. 
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(11) 


Figure  8.  Illustration  of  forces  and  posi¬ 
tion  vectors  used  to  calculate  configurational 
term  in  Irving-Kirkwood  pressure  tensor. 


Figure  9.  Radial  distribution  of  the  normal 
component  of  the  Irving-Kirkwood  pressure 
tensor  for  bare  5nm  gold  nanoparticles  at  var¬ 
ious  temperatures. 


r 

where  7  is  the  surface  tension,  r  is  the  nanoparticle  radius, 
and  A P  is  the  difference  in  pressure  between  the  nanopar¬ 
ticle  and  the  surrounding  medium,  in  this  case  a  vacuum. 
Assuming  a  pressure  of  5900  bar  and  a  radius  of  2.5nm 
the  computed  surface  tension  is  about  738  This  re¬ 
sult  compares  favorably  to  the  published  data  for  bulk  gold 
at  the  melting  temperature  of  742  Another  calcula¬ 
tion  from  Widom  et  al  [22]  is  used  to  calculate  the  surface 
tension  from  the  MD  simulation  with  equation  eq:  widom. 


T2 


Prif 


i  dpN  (r) 
dr 


dr 


(12) 


In  equation  12,  pr2  and  pri  are  the  computed  pressure  at  the 
maximum  and  minimum  model  radii,  respectively.  dp^r^ 
is  the  derivative  of  the  normal  pressure  at  radius  r.  This 
model  has  been  computed  for  the  bare  5nm  gold  nanoparti¬ 
cle  at  900  K  and  1000  K,  and  result  in  surface  tension  values 
of  751  and  617  respectively.  These  temperature 
bound  the  melting  temperature  and  the  surface  tension  re¬ 
sults  bound  the  published  data  of  742  again  providing 
some  validation  of  the  MD  model  used. 


4.  Conclusions 

This  overview  of  simulation  results  for  the  alkanethio- 
late  SAM  coated  gold  nanoparticle  has  provided  some  in¬ 
sights  into  this  system.  For  instance,  from  the  diffusivity 
calculations,  three  distinct  temperature  ranges  have  been 
observed.  These  ranges  could  be  different  phases  for  the 
adsorbed  SAMs,  much  like  those  observed  on  flat  surfaces 
[23],  The  radial  density  results  have  shown  that  the  surface 
of  the  gold  nanoparticle  becomes  corrugated  as  SAMs  are 
adsorbed  onto  the  surface,  and  a  corrugation  factor  has  been 
defined  herein.  These  results  have  spurred  closer  inspection 
of  the  material  system  so  that  more  properties  and  a  phase 
diagram  can  be  developed. 
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